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" Long-wave  infrared  semiconductor  negative  refraction 
metamaterials  for  high-resolution  imaging" 


Executive  summary: 

Goals  of  the  research: 

Metamaterials  form  a  new  paradigm  for  design  and  development  of  novel  optical  and  optoelectronics 
components[l-6].  A  new  type  of  metamaterials,  semiconductor  metamaterials,  operating  at  long-wave- 
IR  frequency  range,  has  been  recently  proposed  and  experimentally  verified  by  the  team  [6-8].  In 
comparison  with  other  metamaterial  designs,  semiconductor  structures  offer  the  advantages  of 

(i)  straightforward  fabrication  that  does  not  require  three-dimensional  patterning 

(ii)  material  compatibility  with  semiconductor  technology, 

(iii)  relatively  low  loss  that  can  be  further  mitigated  by  electric  pumping  via  quantum  cascade 
(QC)  technology 

The  goals  of  the  project  were  to  assess  the  perspectives  of  semiconductor  metamaterials  for  broadband 
high  resolution  imaging.  Specifically,  the  project  aimed  to 

(1)  increase  operating  bandwidth  of  the  metamaterials  by  combining  several  layered  structures 
with  carefully  designed  dispersion 

(2)  reduce  the  loss  of  semiconductor  metamaterials  by  incorporating  quantum  cascade  gain  regions 
into  the  structure,  and 

(3)  develop  a  planar  focusing  system  with  sub-wavelength  focal  spot 

Main  results: 

With  relation  to  the  tasks  above,  we  have 

(1)  theoretically  developed  and  experimentally  verified  a  set  of  multilayered  metamaterials  with 
negative  refraction  response  spanning  the  frequency  range  8... 11  fim  (frequency-independent 
index  spanning  8. ..9. 5  /im),  leading  to  the  increase  of  bandwidth  of  negative  refraction  by  27% 

(2)  designed  and  fabricated  quantum  cascade  structures  integrated  into  semiconductor 
metamaterials;  experimental  characterization  indicates  QC  action  spectrally  overlapping  with 
negative  refraction  region 

(3)  theoretically  developed  a  system  that  is  capable  of  achieving  focal  spots  of  the  order  of  1  jum, 
and  a  new  analytical  technique  that  can  be  utilized  to  image  subwavelength  objects  in  the  far 
field. 
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Goal  1:  enhancement  of  bandwidth  of  negative  refraction  of  metamaterial 

The  basic  principle  of  operation  of  our  semiconductor  metamaterial  lies  in  the  collective  response  of 
highly-doped  plasmonic  layers,  interlaced  by  undoped  layers,  to  the  incident  long-wave  infrared  (LWIR) 
light.  Since  the  thickness  of  individual  layer  (~80  nm)  in  this  structure  is  much  smaller  than  the 
wavelength  of  light  (~8  uni),  the  optical  properties  of  the  composite  can  be  described  by  the  effective 
medium  theory.  In  fact,  the  collective  response  of  the  multilayer  is  identical  to  the  response  of  a 
homogeneous  slab  with  thickness  equal  to  the  total  thickness  of  the  multilayer  with  anisotropic 
dielectric  permittivity  with  components  of  permittivity  tensor  given  by: 

a 1e1+aieI . 

11  a±+a2 

_  CLld2ele2 
1  a1e2+a2e1 

where  a1(  a2  are  the  thicknesses  of  the  components  of  the  metamaterial,  e1,  e2  are  their  permittivities, 
and  £||  and  represent  the  effective  permittivity  of  the  metamaterial  along  and  perpendicular  to  the 
surface  of  the  layers,  respectively. 

The  typical  dispersion  of  metamaterial  is  shown  in  Fig.l. 

Note  that  the  negative  refraction,  signified  by  the 
condition  6||  >  0, 6±  <  0,  is  limited  to  the  frequencies 
between  the  frequency  where  e1  —  0  and  where  £||  = 

0 .  Also  note  that  the  frequency  response  of  the 
metamaterial  is  highly  dispersive,  i.e.  effective 
permittivities  of  the  metamaterial  strongly  depend  on 
frequency.  Although  the  design  of  the  metamaterial 
allows  the  control  of  dispersion  by  adjusting  the  doping 
of  plasmonic  layers,  and  thus  shifting  the  negative 
refraction  band,  the  response  of  the  two-component 
system  remains  dispersive. 

In  order  to  minimize  frequency  dispersion  and  to 
enhance  the  negative  refraction  band  we  designed, 
fabricated,  and  verified  a  multi-component 
metamaterial,  schematically  shown  in  Fig.  2.  The  basic 
principle  of  operation  of  this  structure  is  similar  to  that 
of  multi-lens  systems  corrected  for  chromatic  aberration.  The  main  difference,  however,  is  the  ability  of 
the  relatively  thin  metamaterial  "correcting"  layer  to  significantly  adjust  the  dispersion  of  the  much 
thicker  main  layer; 

Two-  and  four-layer  devices  have  been  fabricated  and  their  dispersion  has  been  analyzed  using 
reflection/transmission  measurements,  as  described  in  [6-8].  The  dispersion  of  the  system  is  most 
clearly  visualized  by  the  angular  position  of  the  Brewster's  angle  in  the  negative  refraction  regime;  the 
spectral  width  of  the  this  regime  can  be  assessed  by  measuring  the  spectral  position  of  the  discontinuity 
of  the  Brewster's  angle  (beginning  of  the  negative  refraction  regime)  and  by  abrupt  increase  of 
reflection  of  TE-polarized  waves  (end  of  the  negative-refraction  regime)  [6,9,10]. 

The  spectral  response  of  the  baseline  sample  is  compared  to  that  of  the  two-  and  four-layer  structures 
in  Fig.3.  Note  that  the  spectrally-corrected  structures  provide  increased  bandwidth  and  significantly 
flatter  response  as  compared  to  the  single-layer  system.  Also  note  that  the  response  of  the  final 


Fig-  1:  Dielectric  function,  and  e h  of  a  n+-i- 
n+  heterostructure  calculated  using  the 
effective  medium  theory  and  nd  =  7.5  x  1018 
cm’3.  The  inset  shows  the  layered  structure 
and  relative  orientation  of  the  dielectric 
function  (eL  and  £||). 
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fabricated  samples  closely  matches  the  desired  response  of  theoretically-optimized  structure. 
Quantitatively,  operating  bandwidth  of  the  sample  with  n  metamaterial  (multi)layers,  defined  as 


with  A1  and  A2  being  the  limits  of  the  negative  refraction  region  is,  yields  the  following  results: 

B1  =  26.2% 

B2  =  32.0% 
fi4  =  33.2% 

The  frequency-corrective  coating  thus  provides  the  increase  of  operating  bandwidth  is  27%;  with  flat 
negative  refraction  regime  observed  between  8  and  9.5  microns  (full  negative  refraction  region 
spanning  8...11 /zm  range)  [10].  We  anticipate  that  with  additional  optimization  we  will  be  able  to 
further  extend  the  operating  range  and  compensate  the  material  dispersion. 


A,^m 


Fig.2  (a)  schematic  of  the  multilayerwith  dispersion  correction;  the  thin  top  layers  provide  compensation  for  frequency 
dispersion  of  the  thick  bottom  layer;  (b)  dispersion  of  each  individual  layer  in  (a);  (c)  beam  shift  due  to  thick  layer  in  (a) 
alone  (dashed  line)  and  the  combined  beam  shift  in  the  dispersion-compensated  system 
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Fig.3.  Analytically  calculated  (left)  and  experimentally  measured  (right)  dispersion  of  single-stack  MM  (top),  two-stack  system 
(middle)  and  four-stack  metamaterial  (bottom).  Frequency  dispersion  is  visualized  by  frequency/angular  dependence  of  the 
Brewster's  angle,  corresponding  to  the  minimum  in  TM-polarized  reflection.  Negative  refraction  region  starts  from  discontinuity 
of  the  Brewster  angle  (~8 fim)  and  ends  when  Brewster  angle  disappears  ( — 11  fim). 
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Goal  2:  loss  reduction  by  incorporating  the  material  gain 


As  mentioned  above,  the  design  of  the 
metamaterial  represents  a  multilayered 
system  of  plasmonic  (conducting)  and 
dielectric  layers.  The  figure  of  merit  (FOM) 
of  the  original  metamaterial,  defined  as  the 
ratio  of  the  real  part  of  refractive  index  to 
the  imaginary  part,  varied  from  2  to  30 
depending  on  the  angle  of  incidence  and 
wavelength  (Fig.4). 

The  idea  of  the  current  thrust  of  the 
program  was  to  replace  the  passive, 
homogeneous  layers  of  the  system  with 
active  quantum  cascade  multistacks, 
providing  electrically-driven  optical  gain  to 
the  device,  thereby  compensating  the 
losses  in  the  system  and  increasing  FOM,  as 
indicated  by  the  changes  in  transient 
transmittance  (TT),  shown  in  Fig.5. 


1000 

100 


1 
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Fig  4.  Decay  constant  a  and  Figure  of  Merit  (FOM)  of  a  passive 
semiconductor  metamaterial  as  a  function  of  angle  of  incidence 
and  frequency 
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Fig.5  Theoretical  calculation  of  transient  transmission  as  a  function  of  wavelength  and  incident  angle  at  10%  loss  compensation 
[corresponding  to  ~100  cm-1  at  20°] 

In  the  course  of  the  project  we  have  fabricated  several  samples  with  QC  active  regions  replacing  the 
passive  dielectric  structures;  our  first  tests  of  the  system  suggest  that 

•  The  presence  of  QC  regions  does  not  destroy  negative  refraction  response  of  the  system 

•  The  VI  characteristics  of  the  device  are  identical  to  those  of  working  QC  laser  (Fig. 6) 

We  are  in  the  process  of  analyzing  the  data  for  second  generation  of  samples,  where  QC  band  is 
spectrally  aligned  with  strong  anisotropy  to  determine  the  effect  of  the  optical  gain  on  the  FOM. 
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Fig-6  Active  metamaterials,  (a)  schematic  design  of  active  metamaterial  where  intrinsic  layers  are  replaced  with  QCL  structures, 
(b)  IV  characteristics  of  the  test  sample.  The  properties  are  identical  to  "conventional"  QCLs;  (c-d)  Transmission  (c)  and 
reflection  (d)  spectra  of  active  metamaterials  [with  no  pumping]  are  similar  to  those  of  passive  structures 
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Goal  3:  development  of  planar  focusing  system  with  subwavelength  focal 
spots  and  of  a  technique  to  measure  the  subwavelength  objects  in  the  far 
field 

One  of  the  great  benefits  of  strongly  anisotropic  metamaterials  lies  in  their  ability  to  propagate 
subwavelength  information  about  the  objects  into  the  far-field.  This  fact  is  most  clearly  seen  from  the 
dispersion  relation  for  the  plane  wave  that  propagates  inside  anisotropic  metamaterial: 


where  kL  and  k n  represent  components  of  the  wavevector  perpendicular  to  and  parallel  to  the  layers. 
As  seen  from  the  above  equation,  plane  waves  carrying  subwavelength  information,  corresponding  to 
|/C|J  »  co/c  —  2n/X0,  have  real  values  of  k±,  and  thus  propagate  in  the  anisotropic  material.  This 
behavior  is  in  strong  contrast  to  light  propagation  in  conventional  (isotropic)  media,  where  the  waves 
with  large  transverse  components  of  wavevectors  must  have  imaginary  component  of  the  wavevector 
along  the  direction  of  propagation,  resulting  in  the  exponential  decay  of  these  modes  (often  called 
evanescent  modes). 


In  this  work  we  utilize  strong  anisotropy  of  the  metamaterial  to  develop  a  planar  structure  for 
subwavelength  focusing  of  light.  The  idea  is  to  deposit  a  pre-designed  diffraction  grating  on  top  of  the 
strongly  anisotropic  medium  that  would  act  similar  to  the  well-known  Fresnel  lens[ll]  by  shaping  the 
intensity  profile  of  incident  light  and  assuring  the  constructive  interference  of  the  portions  of  incident 
beam  at  the  pre-defined  point  in  the  bulk  of  metamaterial.  Similar  to  Fresnel  lens  approach,  the  grating 
would  break  the  incident  plane  wave  into  a  set  of  Fresnel  zones,  with  the  boundaries  between 
neighboring  Fresnel  zones  corresponding  to  the  phase  shift  of  n  in  the  optical  phase  of  the  signal.  The 
coordinates  of  the  boundary  of  the  m-th  Fresnel  zone  are  given  by  the  relationship: 


e\\f2  +  eLxm  ~  J  e\\f2  +  e±xo  ~  ± 


•  (4) 


where  x0  is  the  coordinate  of  the  first  open  Fresnel  zone,  and  the  sign  at  the  right-hand  side  of  the 
equation  corresponds  to  the  sign  of  e±[12]. 


Fig.7  Subwavelength  focusing  in  hypergratings  based  on  planar  metamaterials[?];  sizes  and  positions  of  the  slits  are  given  by 
Eq.(4);  panels  (a,b)  correspond  to  low-loss  nanowire  (a)  and  nanolayer  (c)  structures;  panels  (c,d)  show  field  profiles  in  focal 
planes;  solid  lines:  in  low-loss  systems,  dashed  lines:  full  loss  structures  (see  [12]  for  more  details). 
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The  focusing  performance  of  such  hypergrating  is  illustrated  in  Fig. 7.  It  is  clearly  seen  that  the  device  is 
capable  of  achieving  sub-micron  focal  spots  in  the  far-field  of  the  grating  at  LWIR  frequency  range. 

Experimental  confirmation  of  the  existence  of  sub-wavelength  focal  spots  is  extremely  challenging  since 
no  LWIR  near-field  microscopy  currently  exists.  In  order  to  get  around  this  fundamental  problem,  we 
have  developed  a  novel  technique  for  far-field  subwavelength  imaging.  Similar  to  hypergrating,  this 
technique  relies  on  a  diffraction  grating,  this  time  deposited  at  the  focal  plane  of  the  device  (Fig. 8).  As 
known  from  standard  diffraction  theory[ll],  the  spectrum  of  the  light  transmitted  through  the  grating 
with  period  A  will  comprise  the  spectrum  of  the  incident  light,  combined  with  the  same  spectrum 
shifted  by  a  multiple  of  Bloch  vector  2n/A.  Note  that  when  the  period  of  the  grating  is  smaller  than  the 
wavelength,  the  Bloch  shift  effectively  converts  evanescent  waves  (|/C|J  >2n/A0)  into  propagating 
waves  (|/C||  —  2nm/A\  <  2n/A0),  which  in  principle  enables  the  detection  of  subwavelength  information 
in  the  far  field. 


The  potential  problem  with  this  technique,  however,  lies  in  the  fact  that  the  transmitted  signal  has  the 
spectra  of  propagating  and  evanescent  radiation  superimposed  on  top  of  each  other.  Similar  problem 
appeared  in  previous  realizations  of  diffraction  grating-based  subwavelength  imaging  [13,14],  where  the 
subwavelength  information  was  limited  to  first  Bloch  zone  (m  =  1),  the  far  field  resolution  was  of  the 
order  of/l0/4,  and  the  distinction  between  the  propagating  m  —  0  and  evanescent  parts  of  the 
spectrum  was  achieved  by  rotating  the  polarization  of  the  incident  light. 


■  -2kA  -kx  kA  2kA  ' 

-3k a  -V&co/C  -{eoi/C  3kA 

kx 


Fig.8  (a)  Schematic  of  the  imaging  configuration;  an  object  is  placed  close  to  the  subwavelength  plasmonic  grating  with  period 
(at  )  and  is  irradiated  by  the  plane  wavewith  incident  angle  ;  the  far-field  intensity  is  measured  inside  the  material  with 
permittivity  ;  (b)  transfer  of  evanescent  information  into  propagating  regime  by  the  m-th  diffraction  order  of  the  grating; 
normal  incidence  is  assumed 


In  the  case  of  hypergrating,  we  need  to  decompose  the  information  from  several  Bloch  orders,  resulting 
in  significant  enhancement  of  far-field  resolution  (~A0/20).  To  accomplish  the  above  goal,  we  vary  the 
direction  (not  polarization)  of  incident  light,  measuring  the  field  profile  at  the  given  location  in  the  far 
field  of  the  device  as  a  function  of  the  incident  angle.  We  then  compute  the  transfer  function  tm  of  the 
diffraction  grating  that  provides  the  ratio  of  the  amplitude  of  the  m-th  Bloch-shifted  component  behind 
the  grating  to  the  original  amplitude  of  the  incident  light.  Finally,  we  use  nonlinear  fit[15]  to  map  the 
analytically  calculated  intensity  of  the  signal 


Hy{X’ZAnc)  = 


TjTjWna(k*,n  ~Ko - kx0 )exp(A/m)rr)exp(//cz {qmn ) z) 


(5) 


to  the  experimentally  measured  data.  This  fitting  procedure  results  in  the  amplitudes  of  the  incident 
signal  a(kx )  that  can  be  used  to  reconstruct  the  image  in  front  of  the  grating. 
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Fig-9  Reconstruction  (blue  solid  curves)  of  single  slits  with  width  (a)  X0/2  ,  (b)  A0/4  ,  (c)  T0/8  ,  and  (d)  /to/20  ;  black  dotted  lines 
represent  the  object  size;  vacuum-to-vacuum  measurement  is  assumed 


To  verify  the  developed  formalism[16],  we  simulated  the  far-field  intensity  distribution,  generated  by 
several  slits  ranging  from  A0/2  to  A0/20,  positioned  behind  the  grating  with  period  A  =  2/l0/3,  used  the 
fitting  technique  to  recover  the  size  of  original  slit  based  on  these  far-field  "measurements".  The  results 
of  this  imaging  procedure  are  illustrated  in  Fig.9. 


Fig.10  reconstruction  of  the  subwavelength  focal  spot  of  the  hypergrating  shown  in  Fig. 7,  based  on  measurements  in  vacuum; 
solid  blue  and  dotted  black  curves  correspond  to  retrieved  and  original  images  respectively 

Similar  technique  can  be  applied  to  recover  the  intensity  distribution  in  the  focal  spot  of  the 
hypergrating.  Corresponding  simulations  are  illustrated  in  Fig.10. 

To  summarize,  we  have 

•  theoretically  developed  and  experimentally  demonstrated  frequency-dispersion  compensation 
in  LWIR  metamaterials  [10], 

•  incorporated  QC  systems  into  metamaterials  with  aim  to  demonstrate  first  electrically-pumped 
active  metamaterial 

•  developed  theoretical  formalism  for  creation  of  sub-wavelength  focal  spots  with  anisotropic 
metamaterials  [12],  and  a  new  technique  for  subwavelength  far-field  imaging  [16,  17] 


Our  research  revealed  a  great  potential  of  semiconductor  metamaterials  for  subwavelength  imaging  and 
focusing. 
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We  analyze  wave  propagation  in  coupled  planar  waveguides,  pointing  specific  attention  to  modal  cross-talk  and 
out-of-plane  scattering  in  quasi-planar  photonics.  An  algorithm  capable  of  accurate  numerical  computation  of 
wave  coupling  in  arrays  of  planar  structures  is  developed  and  illustrated  on  several  examples  of  plasmonic  and 
volumetric  waveguides.  An  analytical  approach  to  reduce  or  completely  eliminate  scattering  and  modal  cross¬ 
talk  in  planar  waveguides  with  anisotropic  materials  is  also  presented.  ©  2009  Optical  Society  of  America 
OCIS  codes:  050.1755,  160.3918,  240.6690,  230.7390. 


1.  INTRODUCTION 

On-chip  communications,  surface  plasmon  optics,  and  Si 
photonics  are  all  examples  of  planar  optics ,  where  the  op¬ 
tical  radiation  is  controllably  guided  on  the  plane  of  a 
photonic  chip.  A  number  of  planar  optical  elements  in¬ 
cluding  lenses,  mirrors,  and  on-chip  waveguides — both 
plasmonic  and  dielectric — have  been  recently  designed, 
fabricated,  and  characterized  [1-12].  However,  with  a  few 
exceptions  [11,13,14],  the  majority  of  recent  studies  fo¬ 
cuses  on  in-plane  propagation  of  light  and  neglects  out-of- 
plane  scattering  of  radiation.  In  this  work  we  analyze  the 
out-of-plane  light  scattering  and  modal  cross-talk  due  to 
effective  index  change  inside  planar  waveguides  and  dem¬ 
onstrate  that  this  scattering  may  substantially  affect  the 
propagation  of  confined  modes  in  complex  planar  systems. 
We  design  a  numerical  approach  to  solve  the  problem  of 
scattering  and  modal  cross-talk  in  planar  or  quasi-planar 
structures  that  contain  anisotropic  elements,  and  present 
the  technique  to  reduce  or  completely  eliminate  the 
scattering  and  cross-talk  with  anisotropic  meta¬ 
materials  [15]. 

The  problem  of  out-of-plane  scattering  has  been  of  con¬ 
sistent  interest  to  the  photonics  community.  Although  a 
number  of  finite-difference  and  finite-elements  tech¬ 
niques,  available  today,  can  successfully  solve  the  prob¬ 
lem  of  scattering  in  relatively  small  geometries  (^SlO  X0, 
with  \0  being  free-space  wavelength),  analysis  of  wave 
propagation  in  an  extended  system  is  beyond  the  capabili¬ 
ties  of  methods  that  rely  on  finite-size  meshing  of  space/ 
time.  One  of  the  ways  to  reduce  memory  requirements  to 
calculate  the  field  in  an  extended  structure  with  a  moder¬ 
ate  number  of  scattering  interfaces  is  to  implement  some 
sort  of  wave-matching  technique  where  the  modal  spec¬ 
trum  is  constructed  to  satisfy  the  solutions  of  Maxwell’s 
equations  in  the  space,  and  only  boundary  conditions 
at  scattering  interfaces  are  enforced,  resulting  in 
calculations  of  amplitudes  of  the  modes.  Effec¬ 
tively,  modal  expansion  can  replace  the  need  to  calculate 
all  field  components  at  every  point  of  space  with  the 


need  to  calculate  modal  amplitudes  in  every  region  of 
space. 

One  of  the  first  descriptions  of  the  wave-matching  ap¬ 
proach  and  its  applications  for  highly  conductive  plas¬ 
monic  guides  can  be  found  in  [13,16].  Reference  [13]  also 
describes  scattering  by  planar  guides  with  highly  sym¬ 
metric  cross-sections.  Green’s  function  formalism  has 
been  utilized  to  analyze  out-of-plane  scattering  of  plas¬ 
monic  guides  in  [17].  An  approach  to  calculate  the  modal 
cross-talk  and  scattering  in  ID  guides  was  developed  in 
[18].  Scattering  by  periodically  corrugated  systems  has 
been  analyzed  with  rigorous  coupled  wave  analysis 
(RCWA)  in  [19].  Recently,  the  generalization  of  field  ex¬ 
pansion  to  calculate  scattering  in  plasmonic  planar 
guides  has  been  presented  in  [14].  However,  while  the 
mode-matching  calculations  were  proven  to  be  highly  ef¬ 
ficient,  this  technique  had  failed  in  the  proximity  to  plas¬ 
mon  resonance  condition,  when  the  field  of  a  surface  wave 
is  highly  confined  to  the  proximity  of  a  metal  interface. 
Here  we  present  a  wave-matching  technique  that  is  ca¬ 
pable  of  solving  for  wave  scattering  in  complex  systems 
formed  by  coupled  planar  waveguides. 

The  rest  of  the  manuscript  is  organized  as  follows.  In 
Section  2,  we  present  the  mode  structure  of  an  arbitrary 
planar  guide  used  in  this  work  (Fig.  1).  Section  3  is  de¬ 
voted  to  the  development  of  a  mode-matching  technique 
in  a  quasi-planar  system  comprising  a  uniform  in 
y -direction  array  of  planar  guides.  The  presented  numeri¬ 
cal  approach  is  illustrated  on  examples  of  light  propaga¬ 
tion  in  several  plasmonic  and  metamaterial  systems  in 
Section  4.  Finally,  Section  5  develops  the  formalism  of 
truly  planar  photonics  where  out-of-plane  scattering  and 
modal  cross-talk  are  not  possible  and  presents  an  ap¬ 
proach  to  design  extremely  low-scattering  plasmonic  cir¬ 
cuits. 

Notations.  In  this  work  we  use  the  following  notations 
for  the  electric  and  magnetic  fields  in  the  system:  the  total 
electric  (. E )  and  magnetic  (H)  fields  are  shown  in  italic  let¬ 
ters;  the  fields  of  modal  components  contributing  to  the 
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ometry  of  the  single  multilayer  stack  is  shown  in  (a);  panel  (b)  explains  the  composition  of  top,  bottom,  and  guided  modes;  field  profiles 
in  the  outside  layers  of  the  multilayer  are  shown;  the  interface  between  two  multilayer  stacks  is  shown  in  (c). 


total  field  are  represented  with  symbols  E  and  H;  the  com¬ 
ponents  of  an  individual  mode  in  a  particular  layer  of  the 
multilayer  stack  are  represented  with  calligraphic  sym¬ 
bols  S  and  TL. 

2.  MODAL  SPECTRUM  OF  PLANAR  GUIDES 

We  start  from  analysis  of  a  modal  spectrum  of  a  planar 
waveguide,  schematically  shown  in  Fig.  1(a).  The  struc¬ 
ture  comprises  a  set  of  N  planar  layers  with  layer  inter¬ 
faces  parallel  to  the  yz  plane,  with  the  jth  layer  occupying 
the  space  between  xJ_1<x<xJ,  and  having  (uniaxial)  di¬ 
electric  permittivity  described  by  a  diagonal  tensor  with 


diagonal  components  Sj={sjX  ,sjz  ,sjz}.  We  assume  that  x0 
=  -o°  and  x/v=°°.  In  this  section  we  consider  the  structure 
that  is  infinitely  extended  in  the  yz  plane. 

The  electromagnetic  fields  in  this  layered  system  can 
be  represented  as  a  set  of  transverse-electric  (TE)  and 
transverse-magnetic  (TM)  polarized  waves  (modes).  Each 
mode  of  the  multilayer  constitutes  a  solution  of  Maxwell’s 
equations  that  is  finite  for  -ooc^Cco  [20].  In  a  homoge¬ 
neous  layered  structure,  the  mode  can  be  parameterized 
by  a  combination  of  (i)  its  polarization  (TE/TM),  (ii)  the  in¬ 
plane  components  of  the  wavevector  (ky,kz),  and  (iii)  a  set 
of  layer-specific  complex  coefficients  {aj}  playing  the  role 
of  amplitudes  of  the  mode  components: 


t{ky,kz)  =  < 


ai^i(kx;1(ky9kz),ky,kz)  +  a1i1(-kx.1(ky9kz),ky9kz),  x<x± 


afij{kx.j(kyX),ky,kz)  +  dj  Sj(-  kx.j(ky,kz),ky,kz),  Xj_!  <x<xj 


a+N£N(kx.N(ky,kz),ky,kz)  +  aNSN(-  kx.N(ky,kz),ky,kz) ,  xN_ 1  <  x 


(1) 


H  (ky,kz)  =  < 


ai'Hi{kx.i(ky,kz),ky,kz)  +  k^iikyjkz) ,ky,kz) ,  x<Xi 


aj  Hj{kx.j{ky  ,  kz)  ,  ky  ,  kz)  +  dj  Hj(-  kx.j(ky,kZ),ky,kZ)  ,  Xj_±  <  X  <  Xj 


aN'H'N(K;N(ky,h),ky,kz)  +  kx.N(ky,kz) ,ky,kz) ,  XN_X  <  X 


(2) 


where  the  fields  and  dispersion  relations  in  each  layer  are  described  by 


g-ia)t+ik-f 

'jky+k* 


TE  -  polarized  waves:  < 


g-iait+ik-fg 


Kj  =  - 


+  kl 


{0,  kz,  -  ky} 

y~  iky  +  kz),  kx.jky,  kx-jkz\ 


k2x.J  +  k2y+k!  co- 2 


(3) 
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TM  -  polarized  waves: 


Hi 


n-io)t+ik-r 


k*  +  ki 


(K  +  K)c 


Kj  K  +  K 

- 1 - 

VZ  XX 

Ej  Ej 


ky+kl 


{0, 


exxk 

bj  Kz 


A 


-e?ky} 


«T  J 


(4) 


Note  that  for  a  given  mode  only  two  of  the  amplitudes  aj 
are  independent  of  each  other.  Indeed,  the  remaining  am¬ 
plitudes  can  be  calculated  using  the  well-known  transfer 
matrix  method  [20]: 

(1+W  d-WlM 

(i  -w  a+wJW/’ 

with  polarization-dependent  parameters  <p,-,  aj,  and  Kj 
given  by  tf=exp[i(kx.j+1±kx.j)xj],  aj  =  1/2,  aj 
=  ef'/(2e^1),  Kf=kx.jlkx.j+1,  K™=kx.j^+1/(kx.j+1^). 

For  multilayer  systems,  layer-specific  transfer  matrices 
can  be  multiplied  together,  yielding  the  transfer-matrix 
relating  the  fields  in  any  two  layers  of  the  multilayer 
stack,  and  thus  solving  the  problem  of  reflection  and 
transmission  of  a  plane  wave  by  the  multilayer  composite. 
The  singular  solution  that  corresponds  to  nonzero  scat¬ 
tered  waves  on  both  sides  of  the  multilayer  structure 
(aj,a^0)  with  zero  incident  fields  (a{=a^=0)  corre¬ 
sponds  to  the  eigen  (guided)  mode  of  the  stack  [20] . 

For  each  polarization,  full  spectrum  of  modes  sup¬ 
ported  by  the  stack  includes  three  groups  of  waves  [see 
Fig.  1(b)] .  The  first  group  contains  a  discrete  spectrum  of 
guided  modes,  exponentially  decaying  into  first  and  last 
layers  (a{=a^=0,  a^  =  l).  Here  we  characterize  these 
waves  by  the  in-plane  components  of  their  wavevectors 
{ky,kz}. 

The  remaining  groups  of  modes  contain  the  continuum 
of  waves,  known  as  open-waveguide  modes  [13,14]  (bulk 
modes).  The  first  of  these  groups  represents  the  modes 
originated  by  a  plane  wave  that  is  incident  on  the  layered 
structure  from  the  top  layer  (a{  =  0,  a^=l),  while  the  sec¬ 
ond  group  represents  the  wave  incident  on  the  structure 
from  the  bottom  layer  (a{=l,  a^=  0).  Here,  the  modes  of 
the  first  group  (“top  modes”)  are  parameterized  by  the 
real-valued  kx;N  in  the  top  {Nth)  layer,  while  the  “bottom 
modes”  are  parameterized  by  the  real-valued  kx-i  in  the 
bottom  (1st)  layer. 

In  the  limit  of  symmetric  distribution  of  permittivity 
e(jc)  =  e(-jc),  the  spectrum  of  top  and  bottom  modes  pro¬ 
posed  here  is  equivalent  to  the  earlier  proposed  [13]  com¬ 
binations  of  “standing  wave”  modes  with  symmetric  and 
antisymmetric  x-profiles.  However,  in  contrast  to  the  lat¬ 
ter,  the  combination  of  top  and  bottom  modes  is  more  eas¬ 
ily  generalizable  to  the  case  of  nonsymmetric  (such  as 
plasmonic)  planar  guides.  Note  that  in  a  majority  of  pre¬ 
vious  studies  of  plasmonic  structures  [14],  bottom  modes 
were  explicitly  omitted.  As  explained  below,  this  omission 
becomes  crucial  in  the  regime  of  strong  surface  plasmon 


polariton  (SPP)  scattering,  e.g.,  in  proximity  to  SPP  reso¬ 
nance  or  in  plasmonic  step  geometry  (Fig.  3). 

Overall,  the  field  inside  the  guiding  structure  can  be 
written  as 

£  =  2  [A<«)+ t{ky,kf)  +  A(«r 

Q 

+  f  [A^\kMky,kz(K;N))  +  Atop~0U%, -  kjkx.N)) 

J  0 

+ Abtm+ou  t(kyMKi))  +Ahtm'(kx)  i{ky,-kz{kx,1))-\dkx, 


H  =  2  [A<«)+ E(ky,k^)  +A^~M(ky-  &<9))] 
q 

+  f  [AtoP+(^)H(^,^(^))+At0P_(^)H(^,-^(^)) 

J  o 

+  Aht™\kx)M(ky,kz(kx;l)) 

+  AUm~(kx)my-kz(kx,1))]dkx.  (6) 

Here  we  assume  that  all  the  modes  in  the  layered  ma¬ 
terial  have  the  same  value  of  ky.  This  assumption  does  not 
limit  the  applicability  of  the  developed  technique,  since 
due  to  translational  symmetry,  any  solution  of  Maxwell’s 
equations  in  the  set  of  coupled  waveguides  can  be  repre¬ 
sented  as  a  linear  combination  of  solutions  corresponding 
to  a  set  of  ky  values.  Likewise,  we  assume  that  excitation 
and  response  of  the  system  are  monochromatic  [E,H 
ocexp(-/cL>£)].  The  linearity  of  Maxwell’s  equations  makes 
it  possible  to  generalize  the  developed  formalism  for  the 
arbitrary  pulse  excitation  by  representing  the  incident  ra¬ 
diation  by  linear  combination  of  monochromatic  waves. 

Note  that  in  the  process  of  calculating  the  waveguide 
modes,  it  may  be  necessary  to  determine  the  proper  sign 
of  the  kx  (or  kz )  component  of  the  wavevector  in  a  particu¬ 
lar  layer.  If  the  component  of  the  wavevector  has  complex 
value,  this  sign  is  determined  from  the  requirement  for 
the  mode  to  be  finite  in  its  domain.  If  the  wavevector  com¬ 
ponent  is  real,  its  sign  should  be  determined  to  impose 
the  propagation  of  energy  in  the  positive  x  (or  z)  direction 
[21]. 

The  set  of  waveguide  modes  defined  above  allows  the 
introduction  of  the  scalar  product 

E.XH5)^.  (7) 

where  the  dagger  (f)  corresponds  to  the  adjoined  field,  i.e., 
field  of  the  mode  propagating  in  the  reversed  z  direction 
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[13,18].  It  can  be  shown  [13,14,18]  that  in  a  given 
multilayer  (i)  all  TM-polarized  waves  are  orthogonal  to  all 
TE-polarized  waves,  (ii)  the  guided  modes  are  orthogonal 
to  each  other,  and  (iii)  the  top  and  bottom  modes  may 
have  some  coupling,  depending  on  s1  and  s^:  if  si  =  ejv>  the 
top  and  bottom  modes  corresponding  to  the  same  value  of 
kx  are  coupled  to  each  other  and  are  orthogonal  to  all 
other  modes;  if  one  of  the  two  materials  is  lossy  (as  is  usu¬ 
ally  the  case  with  plasmonic  structures),  the  top  and  bot¬ 
tom  modes  are,  as  a  rule,  orthogonal  to  each  other. 

Note  that,  similar  to  what  has  been  suggested  in  [14], 
the  scalar  product  can  be  calculated  analytically,  signifi¬ 
cantly  speeding  up  the  calculation. 


3.  MODE  COUPLING  ACROSS  MULTILAYER 
STACKS 

A.  General  Formalism 

We  now  turn  to  the  main  point  of  this  work — discussion  of 
coupling  of  the  modes  at  the  boundary  between  the  two 
multilayer  structures.  For  simplicity,  we  present  results 
for  the  case  when  the  interface  is  located  at  z  =  0  [Fig. 
1(c)].  Generalization  of  the  technique  for  other  locations  of 
the  interface  is  straightforward. 

We  are  solving  the  classical  scattering  problem:  finding 
the  fields  scattered  by  the  interface  provided  that  the  in¬ 
cident  fields  are  known.  The  incident  fields  are  repre¬ 
sented  by  the  modes  propagating  in  the  +z  direction  on 
the  left-hand  side  of  the  interface  (z  <  0)  and  by  the  modes 
travelling  in  the  -z  direction  on  the  right-hand  side  of  the 
interface  (z>0);  the  scattered  fields  are  represented  by 
the  modes  travelling  in  the  -z  direction  on  the  left-hand 
side  of  the  interface  and  by  the  modes  travelling  in  the  +z 
direction  on  the  right-hand  side  of  the  interface.  The 
modal  representation  [Eqs.  (6)]  reduces  the  scattering 
problem  to  an  arithmetic  task  of  finding  the  coefficients 
A~l  and  A^  as  a  function  of  A£  and  A^,  which  can  be  solved 
by  imposing  the  following  set  of  boundary  conditions: 

A  =  A;  HLy=HRy, 

ELy  =  ER;,  HL=HRr.  (8) 

As  can  be  explicitly  verified,  the  remaining  boundary  con¬ 
ditions  follow  from  Eqs.  (8). 

In  the  case  of  normal  incidence  (ky=0),  TM-  and  TE- 
modes  do  not  couple  to  each  other.  Correspondingly,  in 
this  case  the  first  two  boundary  conditions  in  Eqs.  (8)  de¬ 
scribe  the  reflection,  transmission,  and  scattering  of  TM- 
polarized  waves,  while  the  remaining  two  conditions  de¬ 
scribe  the  optical  properties  of  TE-polarized  modes. 

In  order  to  solve  the  scattering  problem,  Eqs.  (8)  need 
to  be  converted  into  the  set  of  coupled  linear  equations  for 
the  amplitudes  of  the  scattered  modes.  To  achieve  this 
goal,  we  substitute  the  modal  expansion  [Eqs.  (6)]  in  Eqs. 
(8),  and  subsequently  multiply  the  resulting  expressions 
by  the  adjoined  fields  of  left-  and  right-hand  side  modes, 
as  illustrated  below. 
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B.  Numerical  Implementation  of  the  Algorithm 

In  numerical  simulations,  it  is  necessary  to  replace  the 
continuous  integration  over  kx  with  finite  sums.  Thus, 
Eqs.  (6)  becomes 

E  =  2  [A(m)+E(m)+  +  A(m)'E(m)'Vm), 

m 

H  =  2  [A(m)+H(m)+  +  A(m)"H(mrVm),  (9) 

m 

where  =  ^(ky ,  ±  k^) ,  and  similar  for  H;  the  summa¬ 
tion  in  Eqs.  (9)  goes  over  all  modes  (guided,  top,  and  bot¬ 
tom),  and  the  weight  factors  w  are  equal  to  1  for  the 
guided  modes,  and  are  determined  by  the  integration 
method  used  for  top  and  bottom  modes  [22].  Note  that  the 
number  of  modes  on  the  left-hand  side  of  the  interface 
does  not  necessarily  equal  to  the  number  of  modes  on  the 
right-hand  side  of  the  interface. 

Equations  (8)  now  become 

nl 

2  [Aim)Xm)++AimrEfr]wLm) 

m=  1 

Nr 

=  X  [Ajf)+E^)++A^rE^r]^m), 

m- 1 

Nl 

X  [Aim)Xm)++Aim)Xm)’]*4m) 

m=  1  y  y 

Nr 

=  X  [A^H^+A^H^'rf, 

m= 1  y  y 

nl 

X  [Aim)+Eim)++^L,rEimr]wim) 

m=  1  y  y 

Nr 

=  X 

m= 1 

nl 

X  [Af)+H^)++ALmrH^r]w;Lm) 

m= 1 

Nr 

=  X  [A^)+H^)+  +  AjrrH^r]»Bm)-  (10) 

m= 1 

To  solve  for  NR+NL  unknown  amplitudes,  we  multiply 
first  two  equations  in  Eqs.  (10)  by  the  fields  of  TM- 
polarized  modes  and  integrate  the  resulting  products  over 
x;  similarly,  we  multiply  the  two  latter  equations  by  the 
fields  of  TE-polarized  modes  and  perform  the  integration. 
Assuming  that  the  index  m  first  spans  the  TE-polarized 
and  then  TM-polarized  waves,  the  procedure  results  in 
the  following  two  sets  of  matrix  equations: 

rc+  mn  a  (m)+  ,  mn  a  (m)~  _  &+  mn  a  (m)+  .  mn  a  (m)~ 

Al  +  \l  Al  -  '-skr  Ar  +  Ar  > 

<x+  mn  a  (m)+  ,  mn  a  (m)~  _  <x+  mn  a  (m)+  ,  (x-  mn  a  (m)~ 
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HL  L  HL  L  HR  R  HR  R  ? 

(X+  mn/i  (m)+  (m)“  _  <x  +  mn  \  (m)+  mn 

^HL  +D£L  -D£i?  A#  +D£i?  , 

(12) 

where  the  summation  over  repeated  index  m  is  assumed 
and  matrix  elements  are  given  by 


^ll  =  [^^R  1£<RR-  fimR  ^Dtze]  1^L~  ^KR  A)<kl) 

^21  =  [&kr  1(£mR~  fiyw  1^kr\  1(£<rr  1(^mL~^mR  1&w) 
J&22  -  l^rnR  1(^mR~^mR  ^SRr]  H^SKR  1£mL-fimR 

(15) 


5 

w  >A/Je 


n<iV^E 

«  >a^e 


n<iV^E 

«  >a^e 


n<iV^E 

.  (13) 

n  >NrE 

Note  that  the  modes  of  z  <  0  region  are  not  necessarily  or¬ 
thogonal  to  the  modes  in  the  z>0  region.  Thus,  the  ma¬ 
trices  ^L’^viL’^HR’  &  hr  may  have  substantial  nondiago¬ 
nal  components  describing  cross-talk  of  the  modes  across 
the  interface. 

In  fact,  the  above  matrices  are  square  and  invertible 
only  when  NL=NR,  in  which  case  one  of  Eqs.  (11)  and  (12) 
can  provide  the  information  required  to  solve  the  scatter¬ 
ing  problem.  However,  even  in  this  case,  inversion  proce¬ 
dure  may  lead  to  significant  numerical  problems  and  is 
undesirable.  When  NLJ=NR  these  matrices  are  rectangu¬ 
lar  and  thus,  even  theoretically,  cannot  be  inverted.  To 
overcome  this  difficulty,  we  reduce  Eqs.  (11)  and  (12)  to 
the  following  set  of  equations  that  represent  the  generali¬ 
zation  of  transfer-matrix  formalism  for  coupled  wave¬ 
guide  structures: 

AR  =^A£  +  ^A£ 

A#  =9t21AL  +  9t22A£ 

Al  =£iiA r  +  &i2Ar 
.  - — _  - — _  ,  (14) 

Al  =£2iAr  +  £22^ 

where  A\mt  =  \A{^, . . .  ,A^jfi})  }  and 


^±mn  (m) 

£ft{L\R}  ~  W{L\R}^ 


,±mn  _  (m) 

1ft{L\R}  ~  W{L\Ry 


rt±mn  (m) 

^H{L\R}  ~  W{L\RY 


~±mn  _  ( m )  , 
^H{L\R}  ~  W{L\RY 


j  ^myni£dx’ 

f  <£h ^dx, 

J  -  00 

j  KkEKdx’ 

f  rf&ffld*. 

J  —00 

J/m y<dx> 
[  E&H Kdx’ 

I  BS.Ef*. 

f 

J  —00 


^11  -  1&HL  1(^HL  ~  &HL  1^Hl\  1(^HR  ~  &HL  1^Hr) 

^12  =  V^HL  1(^HL  ~  &HL  ^£1,]  1&£R  ~  &HL 

<  _  _  _  _  ^  _  _  _  _ 

£21  “  [^£L  1(^£L  ”  &HL  1^HLi  \^HL  1(^HR  "  &HL  1^Hr) 

^£22-[^£l  1(^hl~^hl  ^£l]  1(^nL  1(^hr~^hl  1^zr) 

(16) 

Finally,  we  combine  Eqs.  (14)  and  arrive  to  the  generali¬ 
zation  of  the  scattering-matrix  technique  that  solves  the 
problem  of  interlayer  coupling: 

J  Al  =  (/  -  ^12^21)  1'^ll  A?  +  (/  -  £12^121)  1^12^l22  A L 

\  Ar  =  (I  -  9121^12)  1y^2l£ll  A^  +  (/  -  9^21^12)  1^22 

(17) 

with  I  being  identity  matrices.  Equations  (17)  represent 
the  main  result  of  this  work. 

Before  illustrating  the  accuracy  of  the  developed  ap¬ 
proach,  we  would  like  to  underline  its  main  advantages 
and  disadvantages.  The  approach  is  ideal  for  calculation 
of  light  propagation  in  extended  structures  with  a  rela¬ 
tively  small  number  of  multilayer  segments;  the  increase 
in  the  number  of  segments  results  in  additional  memory 
use  for  each  given  segment  and,  correspondingly,  mini¬ 
mizes  the  advantages  of  a  wave-matching  approach  over 
finite-difference  and  finite-element  schemes. 

The  developed  technique  provides  an  efficient  solution 
of  the  problem  of  coupling  the  multilayers  with  high  index 
contrast  (high  optical  mode  density  difference)  by  imple¬ 
menting  a  multilayer-dependent  number  of  modes.  How¬ 
ever,  our  calculations  show  that  careful  design  of  the  spec¬ 
trum  of  the  open-waveguide  modes  is  necessary  when  the 
coupling  to  and  from  the  highly  confined  modes  is  calcu¬ 
lated. 

For  guided  modes,  the  approach  allows  straightforward 
calculation  of  intermode  cross-talk  by  calculating  the 
mode-specific  pointing-flux  and  multiplying  it  by  the  am¬ 
plitude  of  the  given  mode  squared;  similarly,  the  approach 
allows  for  easy  calculation  of  emission  directionality, 
naturally  separating  the  fields  produced  by  guided  modes 
from  the  fields  of  open-waveguide  modes  and  separating 
the  fields  of  the  incident  waves  from  the  fields  of  the  scat¬ 
tered  waves. 


4.  INTERGUIDE  COUPLING  IN  PLASMONIC 
AND  METAMATERIAL  SYSTEMS 

We  now  illustrate  the  accuracy  of  the  presented  approach 
on  several  examples  of  plasmonic  and  metamaterial  sys¬ 
tems. 
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A.  Light  Emission  and  Scattering  by  Single-Mode 
Waveguide 

We  first  consider  light  coupling  to  and  from  the  wave¬ 
guide.  As  an  example,  we  use  a  600-nm-thick  Si  wave¬ 
guide  surrounded  by  air,  and  calculate  the  coupling  be¬ 
tween  this  system  and  a  homogeneous  dielectric  at  X0 
=  1.5  jmm;  8§i=  12.12  [23]. 

To  analyze  the  accuracy  of  our  technique  we  assume 
that  the  system  is  excited  by  the  TM2  mode  with  ampli¬ 
tude  of  1,  and  study  the  percentage  of  the  reflected  light 
into  TM0,  TM1?  and  TM2  modes  as  a  function  of  dielectric 
permittivity  of  the  homogeneous  dielectric  [see  inset  in 
Fig.  2(a)].  For  comparison,  we  have  calculated  the  same 
parameters  with  commercial  finite-element-method 
(FEM)  software  [24] .  The  perfect  agreement  between  the 
results  of  our  technique  and  FEM  simulations  is  shown  in 
Fig.  2.  As  expected,  even  for  this  relatively  simple  system, 
the  scattering-matrix  approach  uses  orders-of-magnitude 
less  memory  than  the  FEM  model.  More  importantly, 
scattering-matrix  formalism  is  easily  extendable  to  the 
case  of  multiple  guides  or  multiple  interfaces. 

The  field  matching  obtained  with  our  technique  for  sD 
=  12.12  are  illustrated  in  Fig.  2(c).  Panel  (d)  of  the  same 
figure  illustrates  the  matching  obtained  in  coupling  be¬ 
tween  an  air-Si-air  guide  and  an  anisotropic  hyperbolic 
metamaterial  with  e^=  3.6  +  0.05/,  e^  =  - 12.2  +  1.36/ 
[25-36]. 

B.  Light  Scattering  in  Plasmonic  Systems 

To  further  analyze  the  accuracy  and  limitations  of  the  de¬ 
veloped  field-matching  technique,  we  have  calculated  the 
scattering  from  the  plasmonic  analog  of  the  Si  guide  pre¬ 
sented  above:  a  600-nm-thick  plasmonic  gap  waveguide 
operating  at  \0  =  0.6  jmm.  We  assume  that  cladding  of  this 
guide  is  composed  from  two  gold  plates  with  eAu=-8.94 


-2  0  2  -2  0  2 


zF  nm  z.  jim 

Fig.  2.  (Color  online)  (a)  Light  reflection  in  the  planar  air-Si- 
air  waveguide  coupled  to  a  homogeneous  dielectric;  \0=1.5  jmm, 
waveguide  thickness  d  =  0.6  jmm;  the  geometry  and  profiles  of  the 
waveguide  modes  supported  by  the  system  are  shown  in  the  in¬ 
set;  excitation  by  TM2  mode  is  assumed;  the  graph  shows  the 
comparison  between  the  technique  presented  here  (lines)  and 
FEM  simulations  (dots);  R2  =  0  since  the  symmetry  of  TMX  mode 
is  different  from  that  of  TM0  and  TM2  modes,  (b)  Reflection  in 
plasmonic  gap  (Au-Air-Au)  waveguide;  X0=d  =  0.6  jmm;  thickness 
of  lines  illustrates  the  convergence  of  the  computations.  Panels 
(c)  and  (d)  illustrate  the  field  distributions  in  the  system  (a)  with 
sD=  12.12  (c)  and  e^=3.6  +  0.05i,  eg'  =  -12.2+1.36i  (d). 


+  1.32/  [23].  The  agreement  between  the  mode  matching 
and  FEM  calculations  of  modal  reflectivity  in  plasmonic 
gap  guides  is  shown  is  Fig.  2(b). 

We  have  further  simulated  the  propagation  of  plas¬ 
monic  modes  across  surface  structure  defects.  In  particu¬ 
lar,  we  have  used  our  approach  to  simulate  the  SPP 
propagation  in  “plasmonic  step”  geometry  (Fig.  3).  As  ex¬ 
pected,  when  the  incoming  SPP  is  travelling  on  the  upper 
side  of  the  plasmonic  step,  most  of  the  incident  energy  is 
converted  into  free-space  modes.  In  contrast,  when  the  in¬ 
cident  SPP  is  travelling  on  the  lower  side  of  the  step,  the 
majority  of  energy  is  converted  into  the  reflected  SPP 
wave. 

Our  numerical  simulations  demonstrate  that  at 
near-IR  frequencies  (when  \sm\^>Sd\  the  scattered  field 
can  be  successfully  decomposed  into  “top”  open  waveguide 
modes,  as  suggested  in  [14].  However,  in  the  proximity  of 
SPP  resonance  (sm  —  -sd),  the  inclusion  of  the  “bottom” 
modes  is  necessary  to  adequately  describe  the  optical 
properties  of  the  system. 

Note  that  the  results  of  wave-matching  simulations  are 
almost  identical  to  those  obtained  with  FEM.  Once  again, 
we  underline  that  the  wave-matching  technique  allows 
for  calculation  of  light  propagation  in  much  larger  sys¬ 
tems  than  the  FEM  system  does. 

To  assess  the  convergence  of  our  method,  we  have  per¬ 
formed  a  set  of  simulations  for  each  SPP  structure  de¬ 
scribed  above,  varying  the  configuration  of  spectra  of  top 
and  bottom  modes  (for  simplicity,  equidistant  kx  spectra 
between  500  and  3000  modes  were  used).  As  expected,  our 
simulations  showed  that  it  is  necessary  to  design  the 
spectrum  of  bulk  modes  to  adequately  resolve  the  SPP 
propagating  at  the  z  =  0  interface.  Interestingly,  the  “aver¬ 
aged”  parameters  (such  as  amplitudes  of  reflected  guided 
modes)  are  much  more  sensitive  to  spectrum  variations 
than  the  matching  of  the  boundary  conditions  at  z  =  0  in¬ 
terface,  which  is  often  considered  to  be  an  indication  of  ac¬ 
curacy  of  a  numerical  method.  The  typical  inter-set  varia- 


Fig.  3.  (Color  online)  Scattering  of  the  SPP  propagating  across 
the  Au-Air  step;  panels  (a)  and  (b)  show  transmission,  reflection, 
and  scattering  (S=l-R-T)  of  an  SPP  that  is  incident  on  the 
step;  geometry  of  the  structures  is  shown  in  insets;  lines  corre¬ 
spond  to  the  formalism  developed  in  this  work  (line  thickness 
represents  data  variation  due  to  changes  in  spectra  of  bulk 
modes);  dots  represent  results  from  FEM  simulations;  panels  (c) 
and  (d)  illustrate  the  field  distributions  obtained  from  scattering- 
matrix  (c)  and  FEM  (d)  simulations  for  d  =  0.7  jmm. 
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tion  of  reflection,  transmission,  and  scattering  are 
illustrated  in  Figs.  2  and  3.  The  latter  figure  also  shows 
the  agreement  between  the  field  distribution  obtained 
with  our  generalized  scattering-matrix  formalism  and 
with  FEM. 


5.  TRULY  PLANAR  OPTICS 

We  now  turn  to  the  analysis  of  intermode  coupling  and 
out-of-mode  scattering  in  planar  optics.  Planar  optics,  in 
general  (and  SPP  optics  in  particular),  are  fundamentally 
different  from  their  free-space  counterparts.  Thus,  when 
a  plane  wave  is  incident  on  the  planar  interface,  the  scat¬ 
tered  field  can  be  decomposed  into  one  reflected  plane 
wave  and  one  transmitted  plane  wave.  In  contrast  to  this 
behavior,  when  one  guided  mode  is  incident  on  the  planar 
interface  between  two  waveguides,  it  generates  a  con¬ 
tinuum  of  open-waveguide  modes  in  addition  to  the  (sets 
of)  reflected  and  transmitted  guided  waves. 

As  seen  from  Figs.  4  and  5  and  from  [14],  the  typical 
interface  between  waveguide  systems  leads  to  scattering 
of  about  20%  of  incident  radiation.  Every  attempt  to 
change  the  effective  index  of  the  mode  in  planar  struc¬ 
tures  is  necessarily  accompanied  by  modal  cross-talk  or 
by  out-of-plane  scattering  of  radiation. 

It  is  possible,  however,  to  utilize  anisotropic  metamate¬ 
rials  to  completely  eliminate  the  cross-talk,  and  to  map 
the  familiar  laws  of  3D  optics  to  optics  of  planar  guides. 
The  main  idea  of  planar  optics  lies  in  the  ability  to  guide 
light  along  the  planar  optical  circuit  with  no  out-of-plane 
scattering  or  modal  cross-talk.  In  this  section  we  assume 
that  the  layers  on  both  sides  of  the  interface  are  aligned 
with  each  other  ( xLj=xRj ). 

In  order  to  realize  the  efficient  control  over  pulse  propa¬ 
gation  in  the  plane,  the  two  layered  structures  must 


mode  number  mode  number 


Fig.  4.  (Color  online)  (a,c)  An  interface  between  an  air-Si-air 
waveguide  and  isotropic  air-(e  =  6.06)-air  guide  leads  to  substan¬ 
tial  modal  cross-talk,  polarization  mixing,  and  out-of-plane  scat¬ 
tering;  while  the  interface  between  air-Si-air  waveguide  and  its 
anisotropic  truly  planar  optics  analog  allows  for  ideal  mode 
matching  with  light  steering  capabilities  (b,d);  guided  modes  in 
(c,d)  correspond  to  mode  number  >2000;  the  system  is  excited  by 
a  TM2  guided  mode  propagating  at  the  angle  30°  to  the  z  =  0  in¬ 
terface;  the  amplitudes  of  modes  in  panels  (c,d)  are  normalized  to 
the  amplitude  of  the  incident  mode. 


z ,  pm  z,  pm 


Fig.  5.  (Color  online)  Anisotropic  coatings  (b,  d)  can  signifi¬ 
cantly  reduce  (and  almost  eliminate)  the  scattering  losses  in 
plasmonic  circuits;  the  figure  shows  reflection,  transmission,  and 
scattering  in  conventional  plasmonic  circuit  (a)  and  in  a  plas¬ 
monic  system  where  the  space  x>0,  z>0  filled  with  material 
eyz  =  1,  sxx  =  eD  (b);  field  distribution  for  eD  =  4  (c)  clearly  shows 
that  only  a  fraction  of  the  energy  of  incident  SPP  is  transferred 
into  SPP  on  the  right-hand  side  of  the  interface;  on  the  other 
hand,  the  interface  between  the  anisotropic  system  with  syz=  1, 
sxx=4  (d)  allows  for  substantial  modulation  of  SPP  index,  and  re¬ 
sults  in  almost  perfect  SPP-SPP  coupling  across  az  =  0  interface. 


(i)  have  the  same  number  of  guided  modes,  and  (ii)  pro¬ 
vide  the  ability  to  independently  control  the  index  of  the 
mode  [crucial  for  steering  the  light]  and  modal  profile 
[crucial  for  optimizing  the  overlap  integrals  involved  in  <£ 
and  f)  matrices] . 

As  shown  in  [15]  on  the  example  of  surface  waves, 
these  conditions  can  be  satisfied  when 


n2e%  =  eS;,  (18) 

with  n  being  the  constant  number  that  does  not  depend 
on  layer  number  j. 

As  can  be  explicitly  verified,  when  Eq.  (18)  is  satisfied, 
all  <E  and  S)  matrices  become  diagonal.  Thus,  the  inter¬ 
mode  coupling  is  absent  across  the  interface.  The  inter¬ 
face  remains  completely  transparent  to  TE-polarized 
waves,  while  reflection  and  refraction  of  TM-polarized 
modes  are  controlled  by  the  ratio  of  out-of-plane  permit¬ 
tivities.  The  direction  and  amplitudes  of  the  reflected  and 
refracted  modes  are  related  to  the  direction  and  ampli¬ 
tude  of  the  incident  modes  via  the  following  Snell’s  law, 

sin(0;)  =  sin(6>r)  =  n  sin(^),  (19) 


and  Fresnel  coefficients, 

A(m)~  k\ £)  -  A(g)+  2^ 

A(m)+  ~  k(m)  +  k(m)  ’  ~  jW  +  £(m) ' 


(20) 


The  above  equations  represent  generalization  of  the  for¬ 
malism  of  Ref.  [15]  to  multilayered  guides. 

The  concept  of  truly  2D  optics  is  illustrated  in  Fig.  4  on 
example  of  an  air-Si-air  system  coupled  to  a  metamate¬ 
rials  waveguide.  As  expected,  reflection  of  a  single  mode 
in  a  conventional  planar  system  is  accompanied  not  only 


Thongrattanasiri  et  al. 


Vol.  26,  No.  12/December  2009/J.  Opt.  Soc.  Am.  B  B109 


by  significant  radiation  scattering  and  modal  cross-talk 
but  also  by  cross-polarization  coupling.  In  contrast,  for 
metamaterial  guides  the  single-incident  mode  excites  one 
reflected  wave  and  one  transmitted  wave. 

Low -Scattering  Plasmonics.  An  interesting  extension  of 
truly  planar  optics  is  possible  in  plasmonic  systems.  From 
the  fabricational  standpoint,  it  is  highly  desirable  that 
the  plasmonic  circuit  is  fabricated  on  top  of  common  me¬ 
tallic  substrate. 

Fabrication  of  an  extremely  low-scattering  plasmonic 
circuit  is  possible  when  using  anisotropic  dielectrics  de¬ 
posited  on  noble  metals  in  the  limit  of  visible  and  near-IR 
frequencies,  where  the  permittivity  of  metal  is  much 
larger  than  the  permittivity  of  dielectric  (|em| 

Consider  the  situation  where  the  in-plane  (yz)  compo¬ 
nent  of  permittivity  of  a  dielectric  is  kept  constant  across 
a  plasmonic  system,  and  only  the  out-of-plane  ( xx )  compo¬ 
nent  of  permittivity  of  the  superstrate  is  varied.  In  such  a 
system  the  propagation  constant  of  the  plasmonic  mode, 
approximately  given  by  [37] 


ki  +  ki^i 


(21) 


can  be  effectively  controlled  by  changing  the  parameter 

xx 
sd  • 

In  the  same  limit,  the  exponential  decay  of  the  mode 

into  the  dielectric  [E  ,H  ^  k^x)  ,  kj  — (o/cyjsff  /|em|] 

does  not  depend  on  its  propagation  constant.  The  only 
source  of  out-of-plane  scattering  in  such  a  structure  is  re¬ 
lated  to  a  weak  dependence  of  the  field  profile  in  metal 
[^jH^exp on  8^  .  Km  —  to/c \J\ sm I ( 1  —  8^  /|8m|). 

As  seen  from  Fig.  5,  anisotropy  of  dielectric  super¬ 
strate,  achievable,  for  example,  with  electrooptic  effect 
[37],  provides  orders-of-magnitude  suppression  of  out-of- 
plane  scattering  in  plasmonic  systems  with  respect  to  iso¬ 
tropic  counterparts. 


6.  CONCLUSIONS 

To  conclude,  we  have  developed  a  reliable  numerical  tech¬ 
nique  for  calculation  of  light  propagation  in  planar  guides 
and  in  arrays  of  planar  guides.  We  have  illustrated  the 
developed  formalism  on  examples  of  photonic,  metamate¬ 
rial,  and  plasmonic  guides,  and  presented  an  approach  to 
utilize  anisotropic  metamaterials  for  minimization  and 
elimination  of  modal  cross-talk  in  planar  optical  circuits. 

This  work  was  sponsored  by  the  Office  of  Naval  Re¬ 
search  (ONR)  (grant  No.  N00014-07-1-0457),  the  National 
Science  Foundation  (NSF)  (grant  ECCS-0724763),  and 
the  U.  S.  Air  Force  Office  of  Scientific  Research  (AFOSR) 
(grant  FA9550-09-1-0029). 
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Abstract:  Broad-bandwidth,  low-dispersion,  optical  metamaterials  are  desired.  Reflection 
measurements  show  that,  by  using  multiple-metamaterial  semiconductor  stacks  of  varying 
thickness  and  doping,  bandwidth  is  improved  by  47%  over  a  single-stack  mid-infrared 
metamaterial,  and  dispersion  appears  reduced. 
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1.  Introduction 

Optical  metamaterials  exhibiting  negative  refraction  have  drawn  considerable  interest  in  recent  years  due  to  their 
potential  for  sub  wavelength  imaging  [1]  and  optical  cloaking  [2].  More  generally,  the  synthetic  nature  of  these 
materials  provides  the  potential  for  significant  engineering  of  their  optical  response,  allowing  great  creativity  in 
designing  new  optical  materials  for  a  multitude  of  applications.  Recently  demonstrated,  all-semiconductor,  mid- 
infrared  metamaterials  show  much  promise  due  to  their  relative  low  loss,  3-D  nature,  and  straightforward  fabrication 
[3].  However,  the  spectral  bandwidth  of  negative  refraction  for  a  given  material  and  the  large  dispersion  in  this 
negative  refraction  spectral  region  currently  make  these  materials  unsuitable  for  applications  involving  wide-band 
signals,  such  as  imaging.  Here,  we  demonstrate  that  these  materials  can  be  engineered  to  produce  a  broadband,  low- 
dispersion  optical  response,  by  forming  composite  structures  of  multiple-metamaterial  stacks. 


2.  Material  design 

The  mid-infrared  metamaterials  used  in  this  study  are  composed  of  interleaved  80  nm  thick  layers  of  lattice-matched 
Ino.53Gao.47As  and  Alo.48lno.52As.  These  layers  are  grown  in  a  stack  by  molecular  beam  epitaxy  (MBE)  on  a  low- 
doped  InP  substrate.  The  InGaAs  layers  are  highly  doped  with  Si  atoms  to  create  a  plasma  resonance  of  free  carriers, 
while  the  AlInAs  layers  are  left  undoped.  Conventionally,  the  InGaAs  layers  are  doped  at  a  uniform  level,  providing 
a  single-wavelength  plasma  resonance  [3].  In  this  work,  we  use  multiple  doping  levels  in  the  same  structure,  in  order 
to  broaden  the  region  of  negative  refraction  and  to  reduce  dispersion  in  the  angle  of  refraction.  This  is  achieved  by 
growing  one  metamaterial  stack  in  a  conventional  manner,  followed  by  growing  another  stack  with  different  doping 
directly  on  top  in  the  same  growth  chamber,  and  so  forth.  By  accumulating  multiple  resonances  on  top  of  one 
another  in  the  same  structure,  the  optical  response  can  be  tailored.  There  are  two  variables  used  in  designing  the 
behavior  of  each  stack:  doping  and  stack  thickness.  The  doping  level  determines  the  critical  wavelength,  %,  which  is 
the  spectral  location  of  the  transition  from  positive  to  negative  refraction.  Higher  doping  leads  to  a  shorter  critical 
wavelength.  The  thickness  of  the  stack  determines  the  length  of  interaction  with  the  material. 


(a) 


(b) 


Fig.  1.  (a)  Calculation  of  expected  lateral  (in-plane)  shift  for  light  passing  through  a  composite  metamaterial  consisting  of 
one,  two,  and  four  conventional  metamaterial  stacks,  (b)  Measured  ratio  of  transverse  magnetic/trans verse  electric 
reflectance  (RTm/Rte)  off  of  the  surface  of  these  three  different  metamaterials  at  an  incidence  angle  of  46  degrees. 
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We  designed  three  structures  in  this  study.  The  first  is  a  baseline  single-stack  structure  consisting  of  40  periods 
of  interleaved  InGaAs/AlInAs  layers  grown  on  InP,  with  InGaAs  layers  doped  at  1 0.9x1 018  cm"3,  for  a  total 
epitaxial-layer  thickness  of  6.4  pm.  The  second  sample  is  a  two-stack  metamaterial,  where  the  first  stack  grown  on 
the  InP  substrate  consists  of  10  periods  (1.6  pm  thick)  with  InGaAs  doped  at  9.0xl018  cm'3,  and  the  second  stack 
grown  consists  of  50  periods  (8  pm  thick)  with  InGaAs  doped  at  lO.OxlO18  cm"3.  The  third  sample  is  a  four-stack 
metamaterial,  where  the  first  stack  grown  on  InP  contains  6  periods  (0.96  pm  thick)  with  InGaAs  doped  at  7.7x1 018 
cm"3,  the  second  contains  7  periods  (1.12  pm  thick)  with  InGaAs  doped  at  8.2x1 018  cm"3,  the  third  contains  13 
periods  (2.08  pm  thick)  with  InGaAs  doped  at  9.0xl018  cm"3,  and  the  fourth  contains  50  periods  (8  pm  thick)  with 
InGaAs  doped  at  1 0.0x1 018  cm"3.  By  using  a  transfer  matrix  method  to  compute  the  coefficients  of  full 
electromagnetic  waves  propagating  through  each  of  the  stacks,  the  beam  shift  (distance  light  translates  in  the  plane 
of  the  stacks  as  it  passes  through  the  sample)  can  be  calculated  versus  wavelength,  shown  in  Fig.  1(a).  This  reveals 
the  expected  bandwidth  and  flatness  of  the  response  in  the  negative  refraction  (negative  shift)  region  for  each 
structure.  From  this  plot,  one  can  see  that  the  multi-stack  metamaterials  have  broader  and  flatter  negative  refraction. 

3.  Sample  characterization 

Each  of  the  three  grown  wafers  was  characterized  through  measurements  of  reflection  off  of  the  epitaxial-layer.  The 
measurements  were  performed  with  a  Fourier-transform  infrared  spectrometer  and  its  internal  source,  along  with  a 
liquid  nitrogen-cooled  mercury-cadmium-telluride  detector.  Reflection  measurements  were  taken  versus  angle  of 
incidence  of  the  source  incident  on  the  sample,  with  angles  ranging  from  36°  to  80°  in  2°  increments.  A  polarizer 
was  used  to  measure  transverse  magnetic  (TM)  and  transverse  electric  (TE)  reflection  at  each  angle,  all  versus 
wavelength.  Spectra  were  normalized  to  a  background  measurement.  Incoming  TM  polarized  light  experiences 
negative  refraction,  whereas  TE  polarized  light  does  not.  A  clear  discontinuity  in  the  Brewster  angle  of  TM  to  TE 
polarized  reflection  is  indicative  of  this  effect.  Fig.  1(b)  compares  spectra  at  an  incidence  angle  of  46°  for  all  three 
samples.  Furthermore,  in  Fig.  2,  one  can  see  the  TM  reflection  over  TE  reflection  ratio  for  all  angles  versus 
wavelength  for  each  of  the  three  samples:  (a)  the  single-stack  material,  (b)  the  two-stack  material,  and  (c)  the  four- 
stack  material.  As  we  see  in  these  spectra,  where  the  negative  refraction  region  is  within  dashed  lines,  the  multiple- 
stack  metamaterials  have  a  broader  bandwidth  of  negative  refraction.  The  negative  refraction  region  starts  where 
there  is  a  discontinuity  in  the  Brewster  angle  (%~7.6  pm  in  Fig.  2(b))  and  ends  where  the  TM/TE  reflection  rises 
sharply  (Rtm/Rte  =  0.5);  the  single-stack  has  a  bandwidth  of  21%,  the  two-stack  has  a  bandwidth  of  31%,  and  the 
four-stack  has  a  bandwidth  of  30%,  which  is  an  increase  by  47%  in  relative  bandwidth  over  the  single  stack.  In 
addition,  dispersion  appears  to  be  reduced  in  the  negative  refraction  region  for  the  multiple-stack  structures,  where 
the  Brewster  angle  is  varying  less  quickly  with  wavelength  in  Figs.  2(b)  and  2(c)  than  it  is  in  Fig  2(a). 


(a)  (b)  (c) 


Wauelengrth  fjjmj  Wawelengdii  (pm)  Wavelength  {^im* 

Fig.  2.  Ratio  of  TM  over  TE  surface  reflection  (color  scale)  versus  angle  of  incidence  and  wavelength  for  (a)  a 
conventional  single-stack  mid-infrared  metamaterial,  (b)  a  broadband  metamaterial  consisting  of  two  stacks,  and  (c)  a 
broadband  metamaterial  consisting  of  four  stacks.  The  negative  refraction  region  is  demarcated  by  dashed  lines. 

4.  Conclusions 

All-semiconductor  mid-infrared  metamaterials  have  been  demonstrated  with  enhanced  bandwidth  and  reduced 
dispersion.  Bandwidth  was  increased  from  21%  in  a  conventional  reference  structure  to  31%  in  a  composite 
structure  featuring  two  metamaterials  stacked  in  succession  (a  47%  relative  increase).  Furthermore,  dispersion  in  the 
negative  refraction  region  appeared  reduced  in  the  composite  metamaterial.  The  relatively  simple  fabrication  of 
these  materials  allows  much  potential  for  creative  manipulation  of  the  optical  response. 

This  work  is  supported  in  part  by  AFOSR  (grant  UFA  955 0-0 9-1-0029). 
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Abstract:  We  develop  a  technique  for  determining  the  size  of  subwavelength  focal  spots  without 
near-field  microscopy,  based  on  carefully  designed  gratings  that  convert  the  subwavelength 
information  into  propagating  waves,  far-field  measurements,  and  computer  post-processing. 
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Subwavelength  focusing  of  electromagnetic  waves  is  of  interest  in  broad  class  of  applications  spanning  all  parts  of 
electromagnetic  spectra.  Some  applications  of  subwavelength  focusing  include  nanoimaging,  nanolithography, 
compact  photonic  elements,  and  ultra-potent  sensors.  Metamaterials  and  plasmonic  systems  promise  to  conquer  the 
diffraction  limit  and  offer  numerous  solutions  to  the  problem  of  sub  wavelength  imaging  [1].  However,  while  the 
subwavelength  light  management  in  metamaterial  systems  may  be  free  from  diffraction  limit,  finite-wavelength  of 
light  in  vacuum  provides  substantial  challenges  in  experimental  studies  of  nanophotonics.  Although  near-field 
measurements  [2]  provide  a  way  to  analyze  subwavelength  field  distributions,  these  techniques  suffer  from 
nontrivial  convolution  of  the  field  pattern  with  tip  function.  Moreover,  near-field  microscopy  is  generally  not 
available  at  IR  and  THz  parts  of  the  spectrum.  Here  we  present  a  technique  that  can  be  used  to  determine  the 
resolution  of  sub  wavelength  focusing  systems  with  far-field  measurements. 

The  technique  relies  on  carefully  designed  grating  to  convert  the  information  about  subwavelength  features  into 
the  propagating  waves,  similar  to  the  technique  proposed  in  [3].  However,  in  contrast  to  [3],  our  technique  does  not 
rely  on  resonant  enhancement  of  subwavelength  information  with  plasmonic  structures.  We  illustrate  the  developed 
formalism  on  the  example  of  restoration  of  subwavelength  focal  spot  of  the  anisotropy-based  planar  focusing 
structure  -  hypergrating  [4],  and  present  generalization  for  reconstructing  any  unknown  subwavelength  images  [5]. 

The  fundamental  difference  between  the  diffraction-limited  and  subwavelength  images  is  seen  in  the  wavevector 
space.  The  spectrum  of  the  sub  wavelength  focal  spot  is  dominated  by  high- wavenumber  components  that 
exponentially  decay  away  from  the  focal  spot.  The  grating,  that  plays  the  role  of  image-reconstructing  structure, 
located  at  the  image  plane  [Fig.  1(a)],  converts  the  evanescent  waves  into  the  propagating  waves  as  described  [6]  by 

k'  =  k0  +  m/cA,  (1) 


where  k0  =  oo/c  =  2n/A0  is  the  incident  field  wavevector,  m  is  the  diffraction  order,  kA  =  2n/A  is  the  grating 
wavevector  parallel  to  itself,  A  is  the  grating  period,  and  k'  is  the  diffracted  wavevector.  The  general  design  concept 
that  governs  the  design  of  the  grating  is  to  suppress  all  but  m  =  ±  1  modes  and  shift  the  evanescent  signals  into  the 
propagating  regime.  This  one-to-one  wavevector  relation  makes  imaging  retrieval  possible. 
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Fig.  1 .  (a)  schematic  of  the  image  reconstruction  system;  the  grating  is  located  at  the  image  plane,  the  measurement  is  done  at  the  distance  z0  from 
the  grating;  (b)  spectrum  distribution  of  hypergratings  with  positive  refraction  (exy  <  0,£z  >  0)  at  the  focal  point  composed  of  only 
subwavelength  information  (\kxc/a)\  >  f~£~z)‘,  (c)  the  range  of  high  wavevector  waves  are  shifted  into  the  interval  of  low  wavevector  waves 
(\kxc/a)\  <  y/£air )  with  m  =  +1  diffraction  modes;  the  length  of  arrows  corresponds  to  kAc/a); 

The  spectrum  distribution  of  hypergratings  with  positive  refraction  metamaterial  (exy  <  0,£z  >  0)  [4,7]  is 
shown  in  Fig.  1(b).  Note  that  this  kind  of  material  only  allows  the  propagation  of  p-polarized  waves.  In  this  work, 
we  let  z  axis  be  the  propagating  direction  perpendicular  to  the  alignment  of  the  grating  and  %  and  y  axes  be  the 
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direction  parallel  to  the  grating.  With  this  metamaterial,  the  subwavelength  signal  occupies  only  high  wavevector 
information  which  is  immeasurable  in  far-field  vacuum.  However,  with  the  assistance  of  the  grating  located  at  the 
focal  plane,  we  may  convert  the  evanescent  waves  into  the  propagating  waves  which  are  detectable  in  vacuum  [Fig. 
1(c)].  The  far-field  magnetic  field  is  described  as 

fef^u)/c+kA 

Hy{x,z)=  J  d/cxa_1(/cz)exp[t(a_1(/cx)x  +  /cz,_1(A:a;)x)] 

— V  £air^/ c 

+  J  dkxa+1(kx)exp\i(a+1(kx)x  +  kz+1(kx)x)\  (2) 

-Je^u/c-kA 


where  a±1(kx )  =  kx  ±  kA  and  kz>±1(kx)  =  Jeair(a)/c)2  —  a\x  .  a±1(/cx)  are  unknown  transmission  coefficients 

of  diffracted  waves.  Note  that  the  range  of  wavevector  integration  covers  all  propagating  regimes  and  some  of 
evanescent  regimes.  The  transmission  coefficients  are  determined  from  measuring  of  the  field  intensity  at  a  distance 

z0  from  the  grating  and  then  fitting  it  with  the  intensity  function  I(x,z0 )  =  |//y(x,z0)|  .  There  are  many  numerical 
techniques,  such  as  curve  fitting  and  interpolation  method,  which  would  support  this  fitting  case  [8]. 

When  the  amplitudes  of  the  transmitted  waves  are  experimentally  measured,  we  Fourier  transform  Eq.  (2)  and 
shift  the  wavevector  back  by  kA.  Sub  wavelength  image  information  is  retrieved  by  inverse-Fourier  transformation. 
The  agreement  of  our  retrieval  technique  and  finite-element  method  (FEM)  simulation  [9]  is  shown  in  Fig.  2(a). 

We  confirm  the  robustness  of  the  approach  by  reconstructing  the  image  of  single  slits  with  width  a  =  A0/20, 
A0/4,  and  A0/2  located  in  an  isotropic  material.  This  case  is  different  from  the  previous  one  because  the  spectrum 
distribution  covers  not  only  the  evanescent  regimes,  but  also  the  propagating  regimes.  Therefore,  in  order  not  to  mix 
the  propagating  and  evanescent  modes,  one  may  use  the  rigorously  coupled  wave  analysis  (RCWA)  technique  [10] 
to  design  a  grating  which  can  suppress  m  =  0  mode  and  magnify  m  =  ±  1  modes.  Following  the  procedure,  Fig. 
2(b-d)  compares  the  results  from  our  method  and  FEM  simulations. 


Fig.2.  (a)  image  on  the  focal  plane  of  the  hypergratings;  dashed  and  smooth  curves  correspond  to  “numerical  experiments”  and  our  retrieval 
technique,  respectively;  (b-d)  retrieval  of  single  slits  with  width  a  =  A0/20,  Z0/4,  and  A0/2  in  the  isotropic  material,  respectively; 

To  conclude,  we  proposed  an  approach  to  determine  the  size  of  the  (sub wavelength)  focal  spots  with  far-field 
measurements  and  computer  post-processing.  The  approach,  that  utilizes  a  pre-designed  diffraction  grating  located 
at  the  focal  spot,  is  easily  applicable  to  majority  of  sub  wavelength  focusing  systems,  including  metamaterials  and 
plasmonic  structures. 

This  work  has  been  partially  supported  by  ONR  (grant  #N00014-07-l-0457),  NSF  (grant  #ECCS-0724763),  and 
AFOSR  (grant  #FA9550-09- 1-0029). 
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Abstract:  We  develop  a  computational  technique  for  imaging  with  sub  wavelength  resolution 
based  on  far- field  intensity  measurements. 

High-resolution  imaging  is  of  interest  for  a  broad  class  of  applications  spanning  all  parts  of  the  electromagnetic 
spectrum.  Unfortunately,  conventional  far-field  imaging  is  fundamentally  limited  by  the  free -space  wavelength.  The 
diffraction  limit  can  be  halved  with  structured  illumination  microscopy  where  the  spectrum  of  the  incident  light  is 
effectively  doubled  via  interference!!].  Alternatively,  in  the  far- field  superlens,  part  of  the  evanescent  radiation 
emitted  by  an  object  is  resonantly  enhanced  via  surface  plasmon  polaritons,  and  is  subsequently  converted  into 
propagating  waves  with  a  subwavelength  diffraction  grating[2].  Both  techniques  rely  on  multiple  measurements  and 
numerical  reconstruction  algorithms  to  perform  imaging  of  an  unknown  object  and  achieve  ~  A0/4  resolution.  Here 
we  present  an  approach  capable  of  non-resonant  imaging  with  resolution  on  the  order  of  A0/20  with  far- field 
measurements  [3]. 

The  fundamental  difference  between  the  diffraction- 
limited  and  sub  wavelength  images  is  seen  in  wave  vector 
space.  The  spectrum  of  a  sub  wavelength  focal  spot  is 
dominated  by  high-wavenumber  components  that 
exponentially  decay  away  from  the  focal  spot.  The  grating, 
that  plays  the  role  of  image-reconstructing  structure, 
located  at  the  image  plane,  and  translates  the  spectrum  of 
the  source  according  to 

k'  =  k0  +  mkA,  (1) 

where  k0  <  co/c  =  2n/A0  is  the  incident  field  wavevector, 
m  is  the  diffraction  order,  kA  =  2n/ A  is  the  grating 
wavevector,  A  is  the  grating  period,  and  k'  is  the  diffracted 
wavevector.  Note  that  diffraction  gratings  can  convert  the 
originally  evanescent  information  into  propagating  waves 
which  can  be  measured  in  the  far  field.  Provided  that  the 
far- field  measurements  of  the  same  object  are  performed  for  different  values  of  incident  angle  (different  values  of 
fc0),  the  contributions  of  different  diffraction  orders  to  the  final  intensity  distribution  can  be  separated  from  each 
other,  and  the  original  field  distribution  can  be  calculated.  Typical  examples  of  the  restored  intensity  distribution  of 
sub  wavelength  objects  are  shown  in  Fig.  1. 

This  work  has  been  partially  supported  by  ONR  (grant  #N00014-07-l-0457),  NSF  (grant  #ECCS-0724763),  and 
AFOSR  (grant  #FA9550-09!-0029). 
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Fig.l.  solid  lines:  images  of  single-slit  sources  of  different 
sizes,  obtained  with  the  developed  technique 
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We  develop  an  analytical  technique  for  retrieving  the  size  and  shape  of  sub  wavelength  objects  using 
far-held  measurements.  The  approach  relies  on  subwavelength  diffraction  gratings  scattering 
evanescent  information  into  the  far  held  along  with  a  numerical  algorithm  that  is  capable  of 
deconvoluting  this  information  based  on  the  far-held  intensity  measurements.  Several  examples  are 
presented,  demonstrating  resolution  on  the  order  of  \0/20.  The  developed  method  can  be  used  at 
any  frequency  range,  and  may  become  a  practical  alternative  to  scanning  near-held  microscopy. 
©  2010  American  Institute  of  Physics,  [doi:  10. 1063/1 .3487779] 


High-resolution  imaging  is  of  interest  for  the  broad  class 
of  applications  spanning  ah  parts  of  electromagnetic  spec¬ 
trum.  These  applications  include  real-time  biomolecule  im¬ 
aging,  nanolithography,  high  capacity  optical  data  storage 
systems,  and  compact  integrated  optical  telecom  solutions.1-5 
Electromagnetic  metamaterials  and  plasmonic  systems  prom¬ 
ise  to  conquer  the  diffraction  limit  and  offer  numerous  solu¬ 
tions  to  the  inherent  imaging  constraints  caused  by  the  dif¬ 
fraction  limit.6-9  However,  metamaterials  alone  may  not  be 
perfect  solutions  to  the  problem.  Indeed,  while  the  light  man¬ 
agement  inside  metamaterial  systems  may  be  free  from  dif¬ 
fraction,  the  hnite  wavelength  of  light  in  vacuum  provides 
substantial  challenges  in  the  experimental  studies  of  nano¬ 
photonics  in  free  space. 

Although  near-held  scanning  optical  microscopy10,11 
provides  a  way  to  analyze  sub  wavelength  held  distributions, 
this  technique  suffers  from  nontrivial  convolution  of  the  held 
pattern  with  the  tip  function.  It  is  also  implicitly  slow  and 
does  not  support  real-time  imaging  since  the  evanescent  sig¬ 
nal  is  collected  on  a  point-by-point  basis.  Moreover,  near- 
held  microscopy  is  limited  to  the  optical  and  gigahertz  parts 
of  the  spectrum  and  is  not  available  at  mid-  and  far-IR  and 
terahertz  frequencies.  Another  technique,  called  structured  il¬ 
lumination  microscopy  (SIM),12,13  was  shown  to  improve  the 
resolution  of  conventional  optics  by  a  factor  of  two  by  ana¬ 
lyzing  the  light  transmitted  through  a  wavelength- scale 
diffraction  grating,  under  a  set  of  different  illumination  con¬ 
ditions.  More  recently,  an  approach  called  the  far  held  super¬ 
lens  (FSL)  that  combines  the  resonant  held  enhancement  of 
near-held  superlens6,14  with  the  scattering  of  (subwave¬ 
length)  surface  plasmons  into  the  propagating  spectrum  has 
been  presented  in  Refs.  15  and  16.  Here  we  present  a  tech¬ 
nique  that  combines  the  benehts  of  subwavelength  resolution 
of  the  FSL  with  the  broadband  performance  of  the  SIM  and 
can  be  used  to  provide  real-time  imaging  of  unknown  sub¬ 
wavelength  objects  in  the  far-held  with  resolution  on  the 
order  of  \0/20  but  is  free  of  SIM  and  FSL  limitations. 

We  demonstrate  the  developed  formalism  by  recovering 
the  shape  and  size  of  several  single  and  double  slit-sources, 
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and  suggest  an  application  of  the  approach  for  restoration  of 
the  held  at  the  subwavelength  focal  spot  of  the  anisotropy- 
based  planar  focusing  structure — hypergrating.9  Our  analyti¬ 
cal  results  are  compared  with  hnite-element-method  (FEM) 
simulations17  with  good  agreement. 

The  fundamental  difference  between  the  diffraction- 
limited  and  subwavelength  images  is  seen  in  wavevector 
space.  The  spectrum  of  a  sub  wavelength  object  is  dominated 
by  high-wavenumber  components  that  exponentially  decay 
away  from  the  object.  Measurement  of  these  evanescent 
components  is  impossible  in  the  far-held  regime.  However,  a 
sub  wavelength  grating,  located  at  the  object  plane  [Fig. 
1(a)],  can  potentially  help  the  image-reconstructing  proce¬ 
dure  by  converting  evanescent  waves  into  propagating 
waves.18  The  wavevector  of  the  scattered  wave  qm  would  be 
related  to  the  wavevector  kx  of  the  incident  wave  by 


<lm=(kx-kxo)+mkA >  (!) 

where  kx0  is  the  parallel-to-grating  component  of  the 
wavevector  of  the  incident  held,  m  is  the  integer  diffraction 
order,  A  is  a  (subwavelength)  period  of  the  grating,  and  kA 
=  27t/A  [see  Fig.  1].  In  a  sense,  the  subwavelength  grating 
shifts  portions  of  the  evanescent  spectra  into  the  propagating 
regime  [Fig.  1(b)],  where  the  originally  evanescent  compo¬ 
nents  could  be  detected  with  far-held  measurements.  This 
procedure,  however,  leads  to  overlapping  of  the  original  eva¬ 
nescent  signals  with  originally  propagating  waves;  so  that 
extra  measurements  are  required  to  deconvolute  the  different 
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FIG.  1.  (Color  online)  (a)  Schematic  of  the  imaging  configuration;  an  object 
is  placed  close  to  the  subwavelength  plasmonic  grating  with  period  A  (at 
z  =  z0bj)  and  is  irradiated  by  a  plane  wave  with  incident  angle  0inc;  the  far- 
held  intensity  is  measured  inside  the  material  with  permittivity  e.  (b)  Trans¬ 
fer  of  evanescent  information  into  propagating  regime  by  the  mth  diffraction 
order  of  the  grating;  normal  incidence  is  assumed. 
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parts  of  the  spectra  from  the  far-held  measurements.  Both 
SIM  and  FSL  rely  on  the  rotation  of  polarization  of  incident 
light  to  attain  additional  measurements  to  perform  the  above 
deconvolution.  This  procedure,  however,  is  not  desirable  for 
systems  where  the  properties  of  the  sample  strongly  depend 
on  polarization.  Here  we  propose  to  fix  the  polarization  of 
incident  light,  and  instead  rotate  its  direction  (given  by  the 
angle  of  incidence  #inc)  as  shown  in  Fig.  1(a).  The  set  of 
measurements  corresponding  to  different  values  of  6mc  is 
then  postprocessed  with  a  numerical  algorithm  to  restore  the 
unknown  object. 

Computationally,  the  problem  of  imaging  can  be  reduced 
to  the  problem  of  finding  the  unknown  (complex)  amplitudes 
a(kx )  of  the  Fourier-decomposition  of  the  light  scattered  by 
the  object  at  the  z=z0bj  plane  as  follows: 
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where  the  z-component  of  the  wavevector  kz  is  related  to  kx 
through  the  dispersion  relation.18  It  is  straightforward  to 
show  that  light  propagating  behind  the  grating  is  described 
by 

Hy(x,z,  0inc)  =  22  wna(kx  n  -  kx())tm(kx  n 

m  n 


~  kxo)txp(iqmnx)exp[ikz(q,ml)z],  (3) 

where  qmn=(kxn-kM)+mkh,  kz(qmn)  =  'Je(Jic^qln,  and  e 
is  the  permittivity  of  the  (isotropic)  medium  where  the  field 
is  measured.  The  index  m  corresponds  to  the  summation  over 
different  diffraction  orders  while  the  index  n  represents  the 
numerical  integration  over  kx  with  weighting  factors  wn  [Fig. 
1(b)].19  The  efficiency  of  evanescent-to-propagating  light 
conversion  by  a  diffraction  grating  is  given  by  the  optical 
transfer  function  tm  defined  as 


iwi^x,n  ^xo)  ~ 


a(kx ,n  ~  ^xo) 


(4) 


where  a(kxn-kx 0)  and  a{qmn )  are  amplitudes  in  front  of 
U^obj)  and  behind  (z>z0bj)  lhe  grating,  respectively.  The 
transfer  function  is  assumed  to  be  a  known  quantity;  in  our 
work,  tm  is  computed  by  the  rigorously  coupled  wave  analy¬ 
sis  technique.20  Our  calculations  show  that  the  grating  with 
A  =  2X0/3  provides  substantial  coupling  between  propagating 
modes  and  evanescent  spectrum  covered  by  the  first  seven 
diffraction  orders  which  potentially  yields  a  resolution  of 
X0/20;  better  resolution  may  be  achievable  with  further  op¬ 
timization  of  the  grating  parameters. 

To  determine  the  unknown  transmission  coefficients 
a{kxn-kx o)  we  fit  the  experimentally  measured  intensity  dis¬ 
tribution  in  the  far  field,  Ip-  \hp\2,  for  a  set  of  values  of  6inc  to 
Eq.  (3)  using  the  nonlinear  least-squares  fitting  technique 
based  on  the  Gauss-Newton  algorithm21  by  minimizing  the 
difference 


p 

F=  2  Up  - 1 Hy(xp,zp,  0i„c>p)|2]2,  (5) 

P=  1 

where  the  summation  is  extended  over  P  data  points 
{0inCp,Xp,Zp,hp}.  Assuming  that  kxn  is  smaller  than  a  cutoff 
wave  vector  kxmax,  the  unknown  transmission  coefficients 
can  be  represented  by  following  Taylor  series: 


FIG.  2.  (Color  online)  Reconstruction  (solid  curves)  of  single  slits  with 
width  (a)  X0/2,  (b)  X0/4,  (c)  X0/8,  ana  (d)  X0/20;  dotted  lines  represent  the 
object  size;  vacuum-to-vacuum  measurement  is  assumed. 


a(k)  =  a0  +  axkl  +  . . .  +  ajkj  +  . . .  .  (6) 

Finally,  once  the  spectrum  of  the  source  is  determined,  the 
object  is  recovered  using  Eq.  (2). 

We  now  illustrate  the  developed  formalism  by  consider¬ 
ing  the  imaging  of  several  single-slit  sources  with  widths  a 
=  X0/2,  X0/4,  X0/8,  and  X0/20  and  double  slit- slit  sources 
with  two  X0/4  slits  separated  by  X0/5.  We  assume  that  all 
linear  dimensions  are  normalized  to  the  free-space  wave¬ 
length  X0,  and  use  the  following  parameters  for  the  plas- 
monic  grating:  permittivity  e^  =  -100,  thickness  h  =  X0/ 10, 
metallic  filling  factor  /=  0.5,  and  grating  period  A  =  2X0/3. 
To  simulate  experimental  conditions,  we  calculate  the  inten¬ 
sity  behind  the  grating  using  a  commercial  software  that 
implements  FEM  solutions  of  Maxwell’s  equations,17  and 
use  the  calculated  intensity  as  “experimental”  input  field  Ip. 
The  good  agreement  of  our  retrieval  technique  for  single  and 
double  slits  and  their  original  configurations  is  shown  in 
Figs.  2,  3,  and  4(a). 

We  further  apply  the  developed  formalism  to  recover  the 
field  distribution  in  the  focal  plane  of  the  planar  subwave¬ 
length  focusing  device  based  on  strongly  anisotropic  (hyper¬ 
bolic)  metamaterials,  known  as  the  hypergrating.9  A  good 
agreement  between  the  original  and  recovered  fields  is  found 
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FIG.  3.  (Color  online)  Reconstruction  (solid  curves)  of  single  slits  [as  in 
Fig.  2]  with  a  random  2%  noise  included. 
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FIG.  4.  (Color  online)  (a)  Reconstruction  of  two  X0/4  slits  separated  by 
\0/5  (vacuum-to- vacuum  measurement),  (b)  Reconstruction  of  the  sub¬ 
wavelength  focal  spot  of  the  hypergrating  based  on  measurement  in  vacuum; 
solid  and  dotted  curves  correspond  to  retrieved  and  original  images, 
respectively. 


when  using  a  grating  with  parameters  sg=^12  962.6 
+  /3726.0,22  h  =  X0/ 80,  /=  0.5,  and  A  =  \0/(Ve^+  ^8o)>  corre¬ 
sponding  to  Cu  at  \0  =  20  yam.  Note  that  the  particular  hy¬ 
pergrating  system,  having  = -23.25  +  /3. 65  and  sz= 24.84 
+  /0.69,  used  in  our  simulations  suppresses  small  kx  compo¬ 
nents  of  the  field  (see  Ref.  9  for  more  details),  so  that  the 
image  recovery  can  be  achieved  with  a  single  9inc  measure¬ 
ment. 

Finally,  we  analyze  the  tolerance  of  the  developed  for¬ 
malism  to  imperfections,  present  in  any  experiment,  by  add¬ 
ing  a  random  2%  noise  to  “experimental”  field  distributions 
from  single  slits,  and  averaging  the  field  recovered  upon  sev¬ 
eral  noise  realizations.  Our  simulations  demonstrate  that  the 
recovered  procedure  is  relatively  stable  under  small  pertur¬ 
bations  (Fig.  3).  Noise  tolerance  will  be  the  main  limiting 
factor  of  the  proposed  technique;  it  is  natural  to  anticipate 
that  the  signals  corresponding  to  higher  values  of  parameter 
m  in  Eq.  (4)  would  be  more  sensitive  to  imperfections  in 
measurements,  as  well  as  to  imperfections  in  fabrication  of 
diffraction  gratings. 

In  conclusion,  we  have  proposed  a  technique  capable  of 
retrieving  unknown  subwavelength  images  with  far-field 
measurements  by  using  a  diffraction  grating  with  subwave¬ 
length  period  and  computer  postprocessing.  As  compared 
with  previously  introduced  techniques,  the  proposed  ap¬ 
proach  is  implicitly  nonresonant,  and  thus  provides  a  rela¬ 


tively  large  operating  bandwidth.  It  also  enables  subwave¬ 
length  resolution  (<X0/20),  and  opens  the  door  for  real-time 
imaging. 
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